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ABSTRACT 

We describe a new approach to calculating photon trajectories and gravitational lens- 
ing effects in the strong gravitational field of the Kerr black hole. These techniques are 
applied to explore both the imaging and spectral properties of photons that perform 
multiple orbits of the central mass before escaping to infinity. Viewed at large incli- 
nations, these higher order photons contribute ~ 20% of the total luminosity of the 
system for a Schwarzschild hole, whilst for an extreme Kerr black hole this fraction 
rises to ~ 60%. In more realistic models these photons will be re- absorbed by the disc 
at large distances from the hole, but this returning radiation could provide a physical 
mechanism to resolve the discrepancy between the predicted and observed optical/UV 
colours in AGN. Conversely, at low inclinations, higher order images re-intercept the 
disc plane close to the black hole, so need not be absorbed by the disc if this is within 
the plunging region. These photons form a bright ring carrying approximately 10% of 
the total disc luminosity for a Schwarzchild black hole. The spatial separation between 
the inner edge of the disc and the ring is similar to the size of the event horizon. This is 
resolvable for supermassive black holes with proposed X-ray interferometery missions 
such as MAXIM, so has the potential to provide an observational test of strong field 
gravity. 

Key words: accretion, accretion discs; black hole physics; gravitational lensing; line: 
profiles; radiative transfer; relativity; 



1 INTRODUCTION 

Black holes are the ultimate test of strong gravity, space- 
time so warped that not even light can escape. By definition 
they have no emission (apart from Hawking radiation), yet 
their immense gravitational potential energy can be tapped 
by any infalling material. This can power a luminous ac- 
cretion flow where the emission has its origin close to the 
black hole event horizon, as is seen in many objects including 
Active Galactic Nuclei, Galactic black hole binaries, Ultra- 
Luminous X-ray Sources and Gamma Ray Bursts. Photons 
emitted in this region are subject to general relativist ic ef- 
fects such as light-bending, gravitational lensing and red- 
shift, as well as special relativistic eff ects as the emittin g 
material will be moving rapidly (e.g. iFabian et all 1200(1) . 
These are well-understood from a theoretical standpoint, so 
accreting objects provide a natural laboratory to test the 
properties of strong gravitational fields. 

Calculations of the relativistic corrections to photon 
properties have been ongoing for nearly three decades, start- 
ing with the classic work of .Cun ningham ( 1975^) who calcu- 
lated the distortions expected on the spectrum of a geometri- 
cally thin, optically thick, Keplerian accretion disc orbiting 
a Kerr black hole. Interest in these calculations dramati- 
cally increased with the realisation that the accretion disc 



could emit line as well as continuum radiation. Iron Ka flu- 
orescence resulting from X-ray irradiation of the accretion 
disc can give a narrow feature, on which the relativistic dis- 
tortions are much more e asily measured than on the broad 
accretion disc continuum iFabian et al.l (^989|). Since then, 
several groups have developed numerical codes that are ca- 
pable of determinin g these effects both for standard discs 
(|Dovciak. Karas Y aaoob 2004) and alternative emission 
geometries (jBurse^ al.ll2004fl . 

While the problem is well-defined, there are many tech- 
nical and numerical issues which arise in calculating the 
effects of strong gravity. Lightbending can result in lens- 
ing which strongly amplifies the emission, so a very small 
area of the accretion flow can make a large contribution 
on the observed flux. Here we describe how to efficiently 
calculate the effects of strong lightbending, and illustrate 
the effectiveness of this approach by using the code to com- 
pute the most extreme lightbending possible , the higher or- 
der images and spectra of an acc r etion disc (IVierg ut jll993|: 
'Bao. Hadrava OstgaardI Il994l : ICadez. Fanton Calvanil 
.1998>) . The first order image is from photons from the under- 
side of the disc which are bent back into the observers line 
of sight, while second order images are produced by pho- 
tons from the upper side of the disc which complete a full 
orbit around the black hole before reaching the observer. 
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Obviously such paths must cross the equatorial plane, so 
are likely to re-intercept the disc. For an optically thick 
disc then this returning radiation adds to the intrinsic disc 
emission, and can enhance the emissivity at small radii 
for extreme Ker r black holes though it has little effect 
for Schwarzchild (ICunning hl^[l97 3; iLaor. Netzer PiranI 
ll99Gl : lAgol Krolikll200n[r 

Nonetheless, the first order image of the disc can dom- 
inate the flux at high inclinations if the optically thick disc 
has rather limited radial extent. Even if it does not, some 
of the higher order image flux can escape to the observer 
through the inner 'hole' in the disc below the radius of the 
minimum stable orbit, a ssuming that the plunging region 
is optically thin (but see lRevnolds BegelmanlllQQTf) . The 
fraction escaping through this inner 'hole' is rather larger for 
Schwarzchild than for Kerr, as the size of the gap between 
the innermost stable orbit and horizon is larger for the non- 
spinning black hole. Obviously, such paths are exquisitely 
sensitive to the gravitational potential, being close to the 
true photon orbit point which is the (unstable) crossover be- 
tween capture by the black hole, and escape to infinity. This 
makes them potentially an excellent test of strong gravity, 
and they could be observable in nearby luminous AGN with 
micro-arcsec imaging X -ray interferometers such as MAXIM 



micro- arcsec ima gmg a 
(iGendreau et al^M) 



Such instruments could also observe the spin of a nearby 
supermassiv e black hole simply from the size of the di- 
rect image (lFukuell200i iTakahashiEool . assuming that 
the mass and distance are known. A disc down to the last 
stable orbit extends down to 6 Rg in Schwarzchild but 
only 1.23 Rg in extreme Kerr. Lightbending is stronger 
in Kerr than in Schwarzchild, but the apparent size of 
the 'hole' in the disc still changes by a factor of 3. 
This contrasts with the case where the accretion flow 
has emission down below the last stable orbit, where the 
size of the true black hole shadow is rather simil ar for 

l200nl: 



both Schwarzch ild and Kerr (|Falcke. Melia. A 

iTakahashU Observations of the galactic black hole 

binaries support to the idea that there is a disc down to 
the minimum stable orb it in cer tain, fairl y high luminos- 
ity spectral states (.Ebisawa et a LlflQQS: K ubota et al.lll999l: 
iGierlihski Donel2004D . However, at lower luminosities this 
is probably replaced by a more complex accretion flow which 
may have continuous energy release down to the horizon 
(Naravan k Yi'"l995; Agol & Krolik 20001: iKrolik HawlevI 
k20Q2 : Afshordi Paczvh ski 2003). While the stellar rem- 
nant black holes require nano-arcsecond imaging to resolve, 
nearby luminous AGN also should have a 'hole' in the 
disc from the minimum stable orbit which is accessible to 
MAXIM. 

The paper is split into two parts. In Section 2 we in- 
troduce the relevant theory necessary to perform efficient 
calculation of photon properties in strong gravity, including 
a review the properties of the effective potentials governing 
photon motion. In Section 3, we apply these techniques to 
examine the contribution of orbiting photons to the standard 
accretion disc solution through images of the inner region of 
the accretion flow, generation of fluorescent Iron Ka line 
profiles and the relative roles played by the different types 
of photon paths. Readers not interested in technical details 
should go straight to Section 3. 



2 CALCULATING PHOTON PATHS IN 
STRONG GRAVITY 

2.1 Null Geodesic Equations 

Material in an accretion flow near to a black hole emits ra- 
diation, which is received by a (possibly distant) observer. 
These photons enable the observer to form an image of the 
flow, which in turn determines the measured spectral prop- 
erties (iBeckwith Donell2004D . In the geometric optics ap- 
proximation, the photons follow null geodesies, which in the 
Kerr metric (written in Boyer-Lindquist coordinates with 
G = M = c = 1) are described by ( Chandrasekhar 198^): 



r = p '^^jRx,q{r) 



(1) 



P=o = p-W^\M 

: p"^A"^ [2ar + A (p^ - 2r) cosec^6>] 

p' = i = p-^A-^ (S^ - 2arX) 

where as usual a is the angular momentum of the source 
(normalised to mass so that a = 1 corresponds to the 'ex- 
tremal' case) and dotted quantities denote differentiation 
with respect to some affine parameter. We introduce the 
metric functions A = — 2r + a^, p^ = + cos^ ^, 
— fr^ + a^ ) — a^Asin^^ and the effective potentials 
(IViergutjllOoil : 



Rx^^ (r) = [(r' + a') - aA] ' - A - (A - af] 



Oa,? (0) = q-\- a cos^ ( 



A^cot^6 



(2) 



Here are constants of motion which are related to 
the photons covarian t ang:ular and linear momenta (see 
iGhandrasekha^ Il983l : Ide Felice fc~eti 1999) . The Kerr 
metric is both stationary and axisymmetric so the derived 
set of geodesic equations (QJ are independent of both of the 
coordinates i and 0. The properties of a given geodesic path 
(specified by a A,^ pair) are completely determined by the 
first two of these equations, which can be merged together 
and integrated to form a governing equation: 



/r 



dr 



Consider the radial motion of the geodesic, described by the 
left-hand side of the above. We fix the ends of the radial 
trajectories at re and ro, which leads us to the general form 
for radial motion: 

dr 



dO 



r,X,q 



(re,ro) = ± 



rr 



dr 



r , 

2 r dr 



(4) 



Here, we have denoted the radial turning point (either per- 
ihelion or aphelion) of the path motion by rt and s^'^ can 
take the values —1, +1 dependent on the nature of the path. 
From the discussion of Wilkins (1972), there are no bound 
null geodesies that terminate above the horizon, implying 
that at most there can be one radial turning point along 
the geodesic path. There are therefore three types of radial 
motion that we must consider: 

(i) From re to ro with no radial turning point encoun- 



tered, implying either (s^ = +l,s^ 



-1 



-1) if ^ 



-1) if 



< To 



> To. 
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(ii) From Ve inwards to perihelion at then outwards to 
To, implying that (s^ = — 1, = +l) for rt < re, rt < Vo. 

(iii) From re outwards to aphelion at rt, then inwards to 
ro, implying that (sj; = +1, = +l) for n > Ve, rt > Vo. 

For the latitudinal motion we can, in general, have an 
arbitrary number of turning points occur along the path (un- 
like in the radial case) . This requires a more complex expres- 
sion to describe the contribution of the latitudinal motion. 
We introduce the variable m — cosO: 

dm 



m,X,q 



(me, mo) 



dm 



\/MA,,(m) 



+ Sr 



1 2 /-, , 2 



J m 

a / 

J mo 

J m_ 



^Mx,q{m) 
dm 



, ^Mx,q{m) 
dm 



(5) 



yjMx,q{m) 



The effective potential now takes the form M\^q{m) 



\^ — q) —a^m^. In the above, m+,- denote the 



2/2 

q-\-m [a 

locations of the latitud inal turning points determ ined by so- 
lution of Mx,q{m) = (iRauch BlandfordllQoi) . whilst Um 
denotes the number of such turning points through which 
the path passes and sj^^ can take the values —1, +1. Specif- 
ically, the case of is described by two possibilities: 

(i) If m+ — — m-, then sl^ — +1,— 1 dependent on 
whether the path is initially directed towards the 'north' 
or 'south' poles of the hole respectively. 

(ii) If m+ 7^ — m-, then sin is in the direction of the 
closest of m+,- to me 

The sign of is determined by the number of turning 
points, Um through which the geodesic passes, = _i^m ^ 
Geodesic paths are therefore described by: 



m,X,q 



(me, mo) 







(6) 



This can be solved analytically in terms of elliptic func- 
tions, which is much more efficient than numerical solutions. 
iRauch BlandfordI (|l994) tabulate these functions to de- 
termine the observed coordinate latitude, Oobs for photons 
with momenta {X,q) emitted from a given radius and lati- 
tude whi ch arr i ve at t he observers' radial coordinate. A lter- 
natively, lAgoll (Il99/tl : ICadez. Fanton Calvail (Il998h fix 
one end of the radial path at infinity with some inclination 
and the coordinates at which the other end crosses the equa - 
torial plane. An additional method, due to IViergutzl (|l993f) 
fixes the ends of the photon paths at the required coordi- 
nates, and a minimisation technique is applied to determine 
valid (A, ^) pairs for a given number of orbits of the black 
hole. 

Our approach combin es as pects of those d escribed by 
iRauch BlandfordI (Il994ll and lViergut3 (Il99.1) . We invert 
the reformulated governing equation to obtain the observed 
co-ordinate latitude of the geodesic, using the technique de- 
scribed byiRauch & Blandford (1994): 



J m 



dm 



.1 .2 



o ^/Mxjm) 



12 12 2 

J m_ 



m-^r,X,q 

+ dm 



(re, To) 



^yMx,qM 

dm 
^Mx,q{m) 



(7) 



We apply th e prop erties of the effective potentials described 
bv IViergut j (|l993r> to dramatically reduce the scale of the 
calculation by analytically restricting the range of A and q to 
those which can escape to infinity. Then we search this range 
for those paths which contribute to a given image order at 
the required observed inclination. These geodesies are then 
projected to form an image of the system on the observers 
sky, which is the n used to determined the measured flux 
(described in iBeckwith Done 20041 . The code therefore 
allows the fast calculation of geodesies linking any two points 
in the spacetime that make a specified number orbits of the 
black hole. 



2.2 The Zeroes of the Effective Potentials 

To obtain solutions of the reformulated governing equa- 
tion we turn to the tables of elliptic integrals pro- 
vided by iRauch BlandfordI (|l994r) , as modified by I A go J 
(|l997f) . These tables, when combined with appropriate nu- 
merical techniques allow us, in principle, to determine the 
geodesies that link an arbitrary emitter, (re, me) and ob- 
server, {ro,mo). 

In practice, however, this calculation is far from trivial. 
By specifying the locations of the emitter and receiver, we 
have placed definite restrictions on the values of the angu- 
lar momentum parameters, (A, ^) for which geodesic motion 
between these two locations is even possible. The geodesic 
motion is dependent on the square root of the two effective 
potentials, Rx,q{r), Mx,q{m), which requires that these func- 
tions remain positive definite at all points along the path. 
If, at any point on the geodesic, this requirement is broken, 
then a potential b arrier is necessa rily formed and no such 
motion is possible. IViergut'3 (|l993l) has shown that these re- 
quirements can be expressed in terms of the interplay of a set 
of algebraic functions, which we consider further here. Note 
that the application of these functions enables us to provide 
tight limits on the region of parameter space which must be 
considered in the calculation and hence hugely reduce the 
scale of the calculation. 

We begin by introducing: 

r — min (re , ro) ; r = max (re , ro) 
m — min (|me| , \mo\) ; m — max (|me| , \mo\) 

The condition that no potential barrier exists between the 
emitter and observer can be re-expressed mathematically as: 



Rx,q (r ^ r ^ r) ^ 0; Mx,q (m ^ |m| ^ m) ^ (8) 

Since the effective potentials are linear in q, we can express 
these requirements as: 



q ^ qr,a (A) ; q^ qm,a (A) 



Here: 



qr,a (A) 



[r^ + a{a-\)\ 



/ A \ A 2 2 2 

qm,a (A) = -A -am 

1 — m^ 



(9) 



(10) 



Physically, these curves correspond to the locus of points on 
the (A, ^) plane for which the given co-ordinate, (r, m) is a 
zero of the associated effective potential. For the latitudinal 
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X \ \ 

Figure 1. Boundary curves on the (A, g)-plane for an extreme (a = 1) Kerr black hole with To = oo, Oe = 7t/2 and Oo = 30°. The 
graph of qm,a (A) is shown by the red curve, qr,a (A) by the green and qc (vc) , Ac (vc) by the blue. On the left-hand panel, the emitter is 
located at re = 6r^, such that re > r^^^ and hence the valid region is bounded by gr,a (A) from above and qm,a (A) from below, whilst 
the qc (tc) , Ac (tc) describes the line below which a photon directed radially inward is inevitably captured by the hole. In the central 
panel, the emitter is now located at re = Srg, such that r^^^ ^ re ^ ^ph^ • "^^^ parameter space is again bounded by qm,a (A) from 
below, however in this case, the upper limit is given by qc (rc) , Ac (rc) for A ^ Ac (re) and by qr,a (A) otherwise. Finally, we see that for 
re = 2r^ (right-hand panel), these boundaries still apply, even though the graph of qr,a (A) is now concave (true for all r ^ 2r^). 



motion, it is found that there are two cases that we must 
consider (taking m as the pericentre of the latidunal motion 
and similarly m is the apocentre): 



Qra, 



, (A) |A| < a^l-m^^ 



(11) 



The description of the radial motion is more com- 
plex, due to the exis tence of the unstable photon orbits 
IChandrasekhaJ (|l983h . These orbits are described by the ex- 
istence of a further set of zeroes of the radial effective poten- 



tial, R\. 



■■ that is subject to the additional constraint 



drRx,q (r) = 0. These conditions yield a pair of paramet- 
ric equations describing a critical curve on the (A,^) plane, 
which define the apparent angular size of the black hole: 



Ac 



a (tc - 1) 



(vr — — rcAc 



a2 (rc - 1)^ 



[4Ac-rc (vc-lf] 



(12) 



The range of values of rc is given by considering the solutions 
of 4a^ = rc (vc — 3)^ (corresponding to the radii of unstable 
photon orbits, eitl ier direct or retrograde , in the equatorial 
plane). Following IChandrasekhaJ (|l983f ). we denote these 



radii by r^^^ and r^^^ and so r^^^ ^ rc ^ r^^\ 

The second complication to this discussion stems from 
IWilkinsI (|l972f) . Specifically, for any particle, orbits with 
E'^ > (where E'^ denotes particle energy; — 1,0 particle 
mass), are unbound, whilst orbits with E'^ < are bound. 
Since photons are massless, this implies that the only bound 
orbits are those with negative E^ (and so must terminate be- 
hind the horizon) , whilst all others must be unbound (except 
for a set of unstable, circular orbits). We cannot therefore 
(unlike in the previous discussion) consider trajectories with 
r as pericentre and r as apocentre. 

For our purposes, three particular cases describing the 
behaviour of the ^r,a(A) curve are of interest. To illustrate 
these, we consider an extreme Kerr (a = 1) hole, with 



the observer located at radial infinity and an inclination of 
Oo = 30°. We locate the emitter in the equatorial plane and 
begin by considering the case r^^'' ^ r (FigureQ] left panel). 
Here, the region of the angular momentum (A, q) plane for 
which no potential barrier is formed between the emitter and 
observer is bounded from below by ^m,a(A), as described by 
equation 1111 The upper limit of this region is provided by 
^r,a(A) and it can be shown that, independent of the loca- 
tions of the emitter and observer, these curves are concave 
and convex functions of A, respectively. This indicates that 
we can provide limits on the allowed range of A via: 



2ar (1 - m^) ± ^/AVT■ 



A+,_ (m, r, a) — 
^{r' + m'a') (13) 



A-a2(l-m2) 

The apparent angular size of the black hole is defined by 
the parametric curve Qc (vc) , Ac (rc). Hence, photons whose 
angular momentum falls under this curve and are directed 
initially inwards towards the hole from r are inevitably cap- 
tured. For r^^'* ^ r, this critical curve intersects the graph 
of qm,a{X) at some Aj'~ (the exact intersection being deter- 
mined numerically) and so inwardly directed photons with 
A^ ^ A ^ AJ, ^m,a(A) ^ ^ Qc{rc), where rc is determined 
by inversion of A = Ac(rc), can be completely excluded from 
the calculation. 

We now move the location of the emitter inwards such 
that r^pl^ ^ r ^ r^^'* ^ r (Figure Q] centre panel). We 
have that the valid region of the angular momentum plane 
is bounded from below by ^m,a(A), as before. The behaviour 
of the upper bound to the angular momentum plane is, how- 
ever, quite different. In this case, since r^^^ ^ r ^ '^I'h' ^ there 
exists some r* = r, which is associated with a (A*,^*) pair, 
through equations (|12|) . For A ^ A*, we therefore have that 
the upper limit of the angular momentum plane is given by 
the critical curve and hence the lower limit on A is deter- 
mined by AJ. Similarly, above A*, the upper limit on the 
plane is given, as in the previous case, by ^r,a(A) and so 
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Figure 2. Closed geodesic loops on the (A, g)-plane (top row) and their projection onto the observers sky (the plane, bottom row) 
for an extreme (a = 1) Kerr black hole with ro = oo, 6e = 7t/2 and 9o = 30°. Geodesies with N = are shown as red loops, those with 
= 1 as blue loops and N = 2 as green loops. Also shown for references on the (A, g)-plane representation are the parameter space 
boundaries (black dashed lines). From left to right, the each panel has re = 8, 6,4, 3, 2, Irg. Note the change in the relative locations of 
the loops as the emitter moves inwards - for re = 8r^, the loops for each image order are completely detached. However, as the emitter 
moves in to 6r^, the = 0, 1 loops overlap and by the time the emitter reaches 4r^, the = loop overlaps both the = 1 and N = 2 
loop. At 3rg, the A^ = loop again crosses the A^ = 1, 2 loops, however in this case the A^ = 1 and N = 2 loops now touch. Moving the 
emitter inward to 2r^, the N = loop detaches itself and moves inside the N = 1,2 loops, which cross. Moving the emitter still further 
inward, to just above r+, we now see that A^ = 1, 2 loops have detached, with the A^ = 1 loop now occurring between the A^ = 0, 2 loops. 




Figure 3. As in Figure [S] for a Schwarzschild (a = 0) black hole, with (from left to right) re = 8, 6, 4, 3, 2.5, 2rp. The loops display a 
qualitatively similar behaviour to that described in the a = 1 case. We note that the restoration of spherical symmetry to the system 
(due to the absence of rotation) has removed the asymmetry of the loops on both planes and hence changed the points of overlap of the 
loops with respect to the a = 1 case. For re ^ 3r^, there are no photons with N = 2 that propagate to infinity, furthermore for re < 3r^, 
there are no A^ = 1 photons that can propagate to infinity either. More importantly, for re < Srg a distant observer is unable to form a 
complete (that is cf) = ^ 27r) image of the ring even using direct (A^ = 0) photons. 



the upper limit on A is given by A+(m, r, a). Note that, for 
A ^ A* there are no photons that are emitted on radially in- 
bound geodesies that reach the observer, whilst for A > A*, 
only those initially ingoing photons with values of q above 
the critical curve can reach the observer. Finally, we note 
that for the emitter located at r ^ 2rg (Figure fig:2.2.1, 
right-hand panel), the graph of ^r,a(A) is now convex. How- 
ever, the relationship between the various curves remains 
unchanged from that described in the previous case and the 
angular momentum plane remains bound. 

2.3 Structure &z Properties of Geodesies Loops 

We begin by considering solutions to the 'governing' equa- 
tion for a simple system where the emitter is located in the 
equatorial plane of an extreme Kerr (a = 1) hole {Oe = 7v/2) 



and takes the form of an infinitesimal ring located at ra- 
dius Te. For such an infinitesimal ring, the flux is undefined 
(since the ring subtends zero solid angle on the observers 
sky) and so we concentrate our attention on the behaviour 
of the solutions on the (A,^) and {a,p) planes (Figure 
To minimise the impact of gravitational lensing, we locate 
the observer at radial infinity and Oo = 30° (as previously). 

Initially, we locate the ring at re — 8rg (Figure [3 panel 
a) and catalogue the complete set of geodesies of zeroth 
(A^ = 0), first {N = 1) and second {N = 2) order. Each 
of these geodesic orders form a closed loop on both the 
(A, q) and (a, (3) planes. We note that the projections on 
both planes are asymmetric about the line A, a = due to 
the breaking of spherical symmetry by the black hole spin 
and also (for the emitter located at 8rg) that the loops for 
each individual image order are completely detached. We 
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Figure 4. Two emitting rings around an extreme Kerr black hole located at re = 6r^ and 6^ = 149° (solid lines), 6^ = 151° (dashed 
lines) with the colour coding as in Figure|3 The observer is located at radial infinity with polar coordinate 6o = 25° (left panel), 6o = 30° 
(centre panel) and Oq = 35° (right panel). In the left panel, \6o\ < |^e'^|, producing circular images of the emitting rings. In the centre 

panel, we now have that |6>J| < |6>o| < \0^\ and the image of the 6^ ring has now separated such that the zeroth and second order images 
appear from the south pole of the hole, whilst the first order image appears from behind the north pole. Finally, in the right-hand panel, 
we have that \6o\ > We'^l and so images of both 6e'^ follow this behaviour. 



now move the emitter inwards to re = Qvg (Figure [21 panel 
b), which causes the loops associated with the zeroth and 
first order images to overlap. The solutions to the geodesic 
equations are therefore multivalu ed at these points in (A, ^) 
space, breaking the statement bv ICunninghamI (|l975l) that 
there are two geodesies linking a point on the accretion disc 
to an observer for each value of the redshift parameter, 
which (for the Keplerian disc considered by Cunningham), 
corresponds to two geodesies at each valid point in A space. 
If we now move the emitter inwards to re = 4:rg (Figure 
121 panel c), we now see that the loops associated with the 
zeroth order geodesic now overlaps with both the first order 
loop and the second order loop. 

We now move the emitter further inwards, so that 
'^ph^ ^ Te = 3rg < r^^'* (Figure 121 panel d). Here, the zeroth 
order loop still overlaps the first and second order loops, 
however, in this case the first and second order loops now 
touch. Moving the emitter inwards still further to re — 2rg 
(Figure 121 panel e), the zeroth order loop now detaches it- 
self from the first and second order loops and moves inside 
these loops on the (a, f3) plane projection. The apparent 
angular size of the zeroth order loop, as measured by the 
distant observer, is now smaller than that of the first and 
second order loops, which has important consequences for 
the calculation of the emergent flux, as we shall shortly see. 
Finally, we move the emitter down to re — rg = 10~^ (Figure 
121 panel f), the location of r+ in Boyer-Lindquist coordinates 
(ho wever, this is not the case in term s of proper distance, 
see lBardeen. Press Teukolskvll 19721) . All of the loops are 
again detached, however, in comparison to the case where 
Ve = 8r^, the ordering of these loops is now reversed on 
the (a,/3) plane, with the second order now subtending the 
greatest angular size of the observers sky. 

We now replace the central extreme Kerr black hole 
with a Schwarzschild (a = 0) black hole and repeat the pre- 
ceding calculation. From Figure El we see that the behaviour 
of the solutions on the two planes is qualitatively similar to 
that of the extreme Kerr case as we move the emitter from 
8rg down to 4r^. However, in this case we note that the 



loops are now symmetric about the line X,a — due to 
the absence of rotation in the system and that this results 
in the quantitative locations of the overlap to change. Note 
that we can understand the existence of these overlaps by 
considering the meaning of the projection of these loops on 
the (A,^) plane. Recall that, for given (A,^) pair, there ex- 
ists a set of roots (n, m^) of the effective potentials Rx^q (r), 
Mx^q (m). Whether or not the geodesic path passes through 
selected members of these sets of roots depends on the ini- 
tial direction taken by the geodesic, described by the Sr, Sm 
parameters in equation |S1 and hence different geodesic paths 
can be described by a single (A,^) pair. 

Consider now the behaviour of these loops as we move 
the emitter to re = Srg for the Schwarzschild hole (Figure El 
panel d). Here, we find that the zeroth and first order loops 
are now detached, but in this case the second order loop 
does not exist. Moving the emitter in further to 2.5r^ (Fig- 
ureEl panel e) , we find that not only does the first order loop 
disappears, but the zeroth order loop is no longer closed! Fi- 
nally, if we move the emitter to re — 2rg = 10~^ (Figure El 
panel f), then we find that we are only able to image a small 
line segment from the ring. To und erstand this behaviour , 
we turn to the discussion given by IChandrasekhaJ ()l983|) 
regarding the 'cone of avoidance', describing the cone gen- 
erated the photons that pass through the unstable circular 
orbits located at 3r^. We let ip denote the half-angle of the 
cone (directed inward at large distances from the hole) and 
it can then be shown that: 

We therefore see that for r ^ 3r^, ^ 7r/2, which implies 
that below 3r^, the apparent angular size of the black hole 
is greater than that of the distant stars. This shows that the 
black hole obscures part of the region = ^ 27r for a 
distant observer when the emitter is located below r = 3r^, 
which implies that the dynamics of the accretion flow in this 
region are extremely difficult to directly measure. 
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Figure 5. The contribution of orbiting photons (higher order images) to a distant observers image of a geometrically thin, optically 
thick, Keplerian accretion disc around Schwazschild (top row) and extreme Kerr (bottom row) black holes. In both cases the observer is 
located at radial infinity with Oo = 15°, the disc extends from the marginally stable orbit (6rg for Schwarzschild, Irg for extreme Kerr) 
to 20rg and the images are coloured by the associated value of the redshift parameter, g = Eo/Ee- From left to right, the panels show 
the contributions from (a) all image orders (A^ = - • • 2), (b) the direct (A^ = 0) image, (c) the first order (A^ = 1) image and (d) the 
second order (N = 2) image. Note that in the Schwarzschild case whilst the orbiting photons are not obscured from the observer by the 
accretion disc, the total fiux (integrated area) measured by the observer from the photons is negligible. Additionally, for the Kerr case, 
the increased radial extent of the disc completely obscures the contribution of the photons. 




Figure 6. As in Figure 13 showing the contribution of photons with N = • • • 2. From left to right, the panels are for 6o = 75.0° • • • 87.5° 
in steps of 2.5°. As the inclination of the system increases, the contribution from photons emitted from the near side of the accretion disc 
directly to the observer is suppressed, reducing the obscuration of the higher order images by the disc as measured by the observer. The 
total area subtended by the unobscured first order (A^ = 1) image increases independently of this effect, due to lensing of photons that 
are emitted from the underside of the far side of the accretion disc, relative to the observer. In the case of the Schwarzschild hole (top 
row), we note that the (total) area of the second order image remains approximately constant with inclination. For the extreme Kerr 
hole (bottom row), photons forming the second order image are completely absorbed by the disc, removing their contribution from the 
overall image received by the observer. Additionally, in contrast to the Schwarzschild case, the images are asymmetric about the a = 
axis due to the effect of black hole spin, particularly noticeable for the higher order cases. For the most extreme inclination, the total 
area subtended by the first order image is on the same order as that from the direct image. 
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2.4 General Emission Geometries 

In the standard picture of accretion onto a massive compact 
object, the emitting material is located in the equatorial 
plane in what is assumed to be a geometrically thin struc- 
ture. As such, gravitational lensing effects only come into 
play for high inclination systems {Oo > 60°). However, if 
one wishes to consider emission from a non-equatorial ge- 
ometry (the case of a geometrically thick, optically thin ac- 
cretion flow, for example), then gravitational lensing effects 
can have important consequences even for low inclination 
observers. As an example, consider again two infinitesimal 
rings, located at re = 6rp, with polar coordinates 0^ — 149° 
and el = 151°. We a gain locate the observer at radial infin- 
ity and initially take Oo = 25° such that \0o\ < \0t'^\ (Figure 
Hleft panel). The two rings form are mapped continuously 
onto the (q:,/3) plane for each individual image order in a 
similar fashion to those considered previously. We now move 
the latitudinal coordinate of the observer to Oo = 30°, such 
that \0^\ < \0o\ < \0l\ (Figure H centre panel). The im- 
age of the lower ring, 0^ is again mapped continuously onto 
the (q:,/3) plane for each individual image order, as would 
be expected as it's relationship to the latitudinal coordi- 
nate of the observer has remained unchanged. This is not 
true however of the image of the upper ring, 0^, which is 
now mapped discontinuously onto the (a, P) plane, with the 
even image orders (zeroth and second) appearing from the 
southern hemisphere of the hole and the odd image orders 
(first) appearing from the northern hemisphere. It is there- 
fore clear that, in order to generate a complete picture of 
the physical properties of such a system, we must include 
the contribution of the higher orders to the calculation. Fi- 
nally, we again move the observer in the latitudinal direction 
to Oo = 35°, such that \0o\ > \0^'^\ (Figure H right panel). 
Now the images of both rings are mapped discontinuously 
onto the (a, f3) plane, which serves to further emphasise the 
importance of the inclusion of the higher order images in the 
calculation. 

3 RESULTS & DISCUSSION 

3.1 Images of Thin Keplerian Accretion Discs 

The contribution of higher order images to the observed flux 
is dependent both on the location of the observer and the 
angular momentum of the hole itself, together with the as- 
sumed geometry and emissivity of the accretion flow. For 
an optically thick accretion disc then any photons which 
re-intersect the disc after emission will be either absorbed 
(and then re-emitted) or reflected by the material. Figure 
[5| shows the contributions of both the direct {N = 0)and 
higher order (N = 1,2) images of a geometrically thin disc 
extending from rms to 20r^, viewed at Oo = 15° for both 
Schwarzschild and maximal Kerr black holes. The principle 
effect of black hole spin for the accretion disk dynamics is 
to move the location of the marginally- stable orbit, rms and 
hence the location of the inner edge of the accretion disc. 
In the case of the Schwarzschild hole, the inner edge of the 
accretion disc is located at 6r^, above the radius of the un- 
stable photon orbits (Srg) so higher order image photons 
which cross the equatorial plane below 6rg are not absorbed 
by the disc and may be able to freely propagate to the ob- 




Figure 7. Unobscured N = 1 (left-hand panel) and N = 2 (right- 
hand panel) images of a thin, Keplerian accretion disc around an 
extreme Kerr black hole, with the disc extending from r-ms to 
20rp . The observer is located at radial infinity with an inclination 
of 6o = 87.5°. The images are highly asymmetric, indicating the 
strong coupling of these photons to the bad hole spin 



server. This contrasts with the extreme Kerr hole behaviour, 
where the accretion disc extends down to Ir^, completely 
obscuring the allowed radial range of the unstable photon 
orbits (Irg ^ Tc ^ 4r^, though mostly 3rp for the polar or- 
bits which dominate the low inclination image) and hence 
the contribution of the higher order photons is completely 
obscured by the accretion disc. 

Figure El shows the same systems viewed at a range of 
large inclination angles, Oo — 75° —87.5°. As in the previous 
discussion for the Schwarzschild hole, the inner edge of the 
accretion disc is located above the location of the photon 
orbits and hence the majority of the orbiting photons are 
able to propagate freely to the observer. However, photons 
in the first order image of the far side of the disc now have 
paths which pass rather closer to the hole than at lower in- 
clination, so the importance of lensing is increased, strongly 
amplifying this part of the image. Most of these photons 
cross the equatorial plane at r ^ 20r^, so can be seen in 
our simulation, but would be obscured by a more physically 
realistic disc which is not entirely flat and has outer edge 
r ^ 20rg. These photons instead would illuminate a large 
region of the underside of the disc as the direct image of 
the disc, adding to its intrinsic emission. By contrast, the 
area on the sky of the second order image remains approx- 
imately constant with increasing inclination, reflecting the 
sensitivity (i.e. instability) of the two-loop photon orbits. 

The high inclination extreme Kerr images are shown in 
Figure El bottom row. The disc itself blocks all the higher 
orbit images close to the black hole, similar to the Oo — 15° 
case. Part of the first order image where the geodesic crosses 
the equatorial plane at r ^ 20r^ can be seen, and this frac- 
tion increases with increasing inclination of the observer. 
Indeed, for the highest inclination system considered in this 
work, the apparent angular size of the first order image is 
approximately equivalent to that subtended by the direct 
image. By contrast to the Schwarzschild case, we note that 
the images of the extreme Kerr system are strongly asym- 
metric about the horizontal axis due to the effect of the 
black hole spin, with the degree of this asymmetry increas- 
ing with the inclination of the observer. This asymmetry is 
most pronounced for the (unobscured) first and second or- 
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der images generated from the highest inclination system 
(see Figure 0. Here, the effect of the range of radii for or- 
biting photons can be seen most clearly in the second order 
image. The image is offset from zero as the ones on the left 
are the photons which are emitted from the side of the disc 
approaching the observer, so go with the spin of the black 
hole and orbit at r^, while photons on the right are the 
retrograde ones at 4rp. The hole in the inner disc is sym- 
metric at r^, so it (just) obscures all the prograde photons 
and easily obscures the retrograde ones. A small decrease in 
spin (e.g. to a=0.998, the equilibrium spin for thin disc ac- 
cretion) mcreases Vrns to 1.23, while the photon orbits span 
1.07 — 3.998r^, so a small part of the prograde higher or- 
der image can be seen. The fraction of higher order photons 
which can escape thro ugh the 'hole' increases with i ncreas - 
mg spm until a < 0.65 lBardeen. Press Teukolskvl 
where the photon orbits are all within Vms- 

3.2 Spectral Properties of Higher Order Images 

To understand how the astrophysical properties of the ac- 
cretion flow couple to the gravitational field of the black 
hole, we g enerate relativstica l ly sme ared line profiles as de- 
scribed in iBeckwith Doii3 (120041) . For clarity, we briefly 
recap this method here. We consider an intrisically narrow 
emission line with rest energy, Ee^ for which the flux distri- 
bution measured by a distant observer at an energy Eo is 
given by: 

Fo (Eo) = ^ J j 9^^ {re, lie) S {Eo - gE^nt) dad(3 (15) 

Here, s(re,/ie) is the local emissivity of the accretion disc, 
which we take to have the form £(re,/ie) = e(re)/(/ie), 
where fie describes the initial direction of the photon rela- 
tive to the local z-axis of the emitting material. The flux at 
each redshift in each image order is then calculated directly 
from the area subtended on the observers sky, together with 
the intrinsic disc emissivity (both radial and angular). This 
gives the transfer function for monochromatic flux, so the 
observed emission is the convolution of this with the intrin- 
sic disc spectrum. In the case of intrinsically monochromatic 
radiation, e.g. the iron Ka fluorescence line which can be 
produced by X-ray illumination of cool, optically thick gas, 
then this transfer function directly gives the expected line 
profile. Photons produced by this emission process in regions 
close to the black hole enable us to examine the properties 
of the multiple orbit photons considered in the preceding 
section from a spectroscopic perspective. 

To generate the line profiles, we apply a radial emis- 
sivity of the form e{re) oc (co nsistent with gravita- 
tiona l energy release within the disc, IZvcki. Done SmithI 
Il999h and consider three possible angular emissivity laws: (i) 
/(/ie) cx: cons., corresponding to an optically thick disc (red 
lines); (ii) /(/ie) cx: /i^^, correspondin g to an optically thin, 
limb brightened disc (blue lines) (iMatt. Perola Stellal 
Il993h and (iii) /(/Xe) oc (1 + 2.06/ie) corresponding 
to an optically thick, limb darkened disc (green lines) 
(jLaor. Netzer PiranI Il990h . Figure |H1 shows the line pro- 
files generated for the low inclination Schwarzschild disc dis- 
cussed in the preceding section. Whilst in the case of the 
extreme Kerr black hole (Figure |^ bottom row), all orbit- 
ing photons returned to the inner region of the disc, for the 



Schwarzschild hole (Figure |51 top row) these photons cross 
the equatorial plane below r ms and can freely propagat e 
to the observer (although see lRevnolds Begelma'nlll997f) . 
Examination of their spectral properties shows that these 
photons play a limited role in forming the overall spectral 
shape measured by a distant observer, which is completely 
dominated by the contribution from the direct image. 

Figure El shows the line profiles generated for a 
Schwarzschild hole, with the accretion disc now viewed at 
high inclination (ignoring the effect of obscuration). For the 
direct image, limb darkening boosts the effects of gravita- 
tional lensing, enhancing the flux from the far side of the 
hole. This is because these photons are strongly bent, i.e. are 
emitted from a lower inclination angle than that at which 
they are observed, s o a limb darkening law m eans that the 
flux here is higher (jBeckwith Pond 120041) . The doppler 
shifts are rather small for this material, so this lensing en- 
hances the flux in the middle of the line. Since the line pro- 
files are normalised to unity, this means that the blue wing 
is less dominant. 

The first order spectra are shown in Figure El middle 
row. The transfer functions mostly retain the characteristic 
double peaked and skewed shape, and again the principle 
effect of the different angular emissivities is to alter the bal- 
ance between the blue wing and lensed middle of the line. 
However, there is some new behaviour for the limb bright- 
ened emissivity. This has the largest change in emissivity 
with angle, and this combined with the exquisite sensitivity 
of lensed paths means that this picks out only a small area 
on the disc, leading to a discrete feature in the spectrum. 
These profiles also show enhancement of the extreme red 
wing of the line, as the photons which orbit generally are 
emitted from the very innermost radii of the disc. 

The discrete features are completely dominant for all 
emissivities at second order (Figure El bottom row). These 
are images of the top of the disc where the photons have 
orbited the black hole, so the paths are even more sensi- 
tive to small changes than first order. Thus the profiles are 
significantly more complex in structure, being dominanted 
by lensing. There are blue and red features at the extreme 
ends of the line profile which are picking out the maximum 
projected velocity of the innermost radii of the disc. These 
have the standard blue peak enhancement. However, the two 
strong features red ward of this are a pair of lensed features, 
from the near and far side of the disc. 

Figure ITHl shows the line profiles generated for the ex- 
treme Kerr hole at high inclinations, again ignoring the effect 
of obscuration. For the direct image (Figure^] top row), 
the lines exhibit the ch aracteristic triangular shape previ- 
ously reported by (e.g) iLaoJ (|l99lh . with the variation in 
angular emissivity acting to alter the balance between the 
different regions of the line on a ^ 5% level. The lines as- 
sociated with the first order image (Figure ITUl middle row) 
exhibit a marked difference in comparison to those associ- 
ated with the Schwarzschild black hole. In general, they are 
broader that those associated with the direct image and re- 
semble a skewed Gaussian combined with a narrow line (due 
to caustic formation) at g ~ 1.0. Here the principle effect of 
changes in the angular emissivity is to alter the height of 
the blue wing, relative to the rest of the line. Again, the line 
profiles associated with the second order image (Figure Cni 
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Figure 8. Relativistic line profiles generated from images (Figure |3 top row) of a geometrically thin, Keplerian accretion disc around a 



Schwarzschild (a = 0) black hole. The observer is located at ( 



15° and the disc extends from rms to 20r^. All profiles are normalised 



such that they contain one photon and are unsmoothed. From left to right, the panels show the contributions from (a) photons with 
= • • • 2, (b) = 0, (c) = 1, (d) = 2. In all cases a e(re) oc r^^ radial emissivity has been applied to the line profile. We consider 
three different types of angular component (i) f{fie) oc cons, (red lines), (ii) /(/ie) oc fi'^^ (green lines) and (iii) f{fie) oc (1 + 2.06/ie) 
(blue lines). We display the lines generated from the direct image (red line), first order (blue lines), second order (green lines) and the 
combined line profile (black line) using a standard £{re) oc r^^. the principle effect of the higher order images is to boost the overall flux 
in the system by ~ 4%, whilst leaving the overall shape of the line approximately unchanged. 
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Figure 9. The contribution of orbiting photons (higher order images) to the Iron Ka line profiles for a geometrically thin, Keplerian 
accretion disc around a Schwarzschild black hole. In all cases, the inner edge of the accretion disc is located at the innermost stable orbit, 
which occurs at 6rg and the outer edge of the disc is located at 20r^ . The observer is located at radial infinity with inclination (from left 
to right) Oo = 75° • • • 87.5° in 2.5° steps. All profiles are normalised to the apparent area on the observers sky and are unsmoothed. The 
line profiles are generated by photons in (a) the zeroth order (direct) image (top row), (b) the first order image (middle row) and (c) 
the second order image (bottom row). Obscuration by the disc is ignored. In all cases a e(re) oc r^^ radial emissivity has been applied 
to the line profile. We consider three different types of angular component (i) f{fie) oc cons, (red lines), (ii) /{fie) oc /j,'^^ (green lines) 
and (iii) f{n-e) oc (1 + 2.06/ie) (blue lines). Line profiles generated by the zeroth order photons have the standard skewed, double peaked 
structure. Those generated by the first order photons have a similar structure, whilst those from the second order photons are far more 
complex. As the observer moves towards higher inclinations, gravitational lensing effects become apparent in the line profile generated 
from zeroth order photons. There is a corresponding increase in the complexity of the profile generated from the first order photons. 
Note however, that the structure of the profile generated from the second order photons remains approximately constant with increasing 
inclination. 
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Figure 10. As in Figure 1771 for accretion discs around an extreme Kerr black hole. By contrast to the Schwarzschild case, the lines 
extend down to lower energies due to the reduced separation between the inner edge of the disc (at the marginally stable orbit) and the 
black hole event horizon. Additionally, the line due to the first order photons extends to higher energies than that due to the zeroth 
order photons, with the upper edge of the line occuring g \n the most extreme case.. Both the first and second order line profiles 

have significantly increased structure compared to the zeroth order profile. The shape of the zeroth and second order line profiles remain 
approximately constant with increasing inclination. However, the internal structure of the first order profile shows a distinctive change, 
with the fraction of radiation in the blue wing of the profile exhibiting a marked increase compared to that in the red wing as the observer 
moves towards the equatorial plane. 



bottom row) are completely dominated by discrete features, 
as in the Schwarzschild case. 



3.3 Image Luminosities 

To understand the relative roles played by each (unobscured) 
image order, we consider the variation of the luminosity of 
each image as a fraction of the total luminosity of the system 
as a function of inclination, again for both Schwarzschild and 
extreme Kerr black holes (Figure^J). These luminosities are 
generated from the integral in redshift space of the line pro- 
files considered in the preceding section and hence we con- 
sider systems with identical properties to those previously 
discussed. In the case of the Schwarzschild hole, we see that, 
for inclinations < 80° , the first order image can be regarded 
as a ^ 10% correction to the emergent flux from the system. 
For inclinations ^ 80°, this image contains 10 — 20% of the 
emergent flux, i.e. even at these high inclinations, the first 
order image can still be regarded as a correction to the di- 
rect image. For the second order image, we see that it plays 
a ^ 1% role independent of inclination. Third and higher or- 
der images will have correspondingly smaller fluxes, so can 
safely be neglected. 

There is a larger fraction of flux in the higher order im- 
ages for the Kerr black hole. At low inclinations, the first 
and second order images contain 10% and 1% of the total 
flux, respectively. The variation of the luminosity fraction 



for the three image orders displays an approximately power 
law like behaviour for inclinations ^ 60°, where a distinct 
break occurs, due to the appearance of a caustic in the first 
order image, whose luminosity is therefore significantly en- 
hanced. It is at this point that the peak luminosity of the 
second order image occurs, which here is on the level of 
10%, approximately an order of magnitude higher than in 
the Schwarzschild case. Remarkably, for inclinations > 75°, 
the luminosity of the first order image produced by the opti- 
cally thick discs is greater than that produced by the direct 
image. 

However, most of the higher image order flux is ex- 
pected to reintercept the disc plane, and hence be absorbed 
and reradiated. This is especially the case for realistic discs 
around an extreme Kerr black hole, where the whole of the 
equatorial plane is covered by the disc from rg to large radii. 
However, for Schwarzchild, the existance of a central 'hole' 
means that the flux from higher order images can escape. A 
realistic disc around a Scharwzchild black hole when viewed 
face on has 10% of its flux in a higher order image ring 
(dropping to ~ 7% for a more extended disk ranging from 
Tms 400r^). Essentially, a spatial resolution equivalent 
to 2rg is required to resolve these features, which for typ- 
ical nearby Active Galactic Nuclei corresponds to an an- 
gular resolution of 0.01 microarcseconds. This is obviously 
extremely technically difficult, but is feasible for an X-ray 
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Figure 11. Variation of (unobscured, integrated) luminosity as- 
sociated with each image as a fraction of the total luminosity 
of the system with inclination for a thin, Keplerian accretion 
disc around a Schwarzschild black hole (left panel) and extreme 
Kerr (right panel). In both cases, the disc extends from rms to 
20rg and the observer is located at radial infinity. We show the 
contributions from (a) the direct image (red lines), (b) the first 
order image (blue lines) and (c) the second order image (green 
lines). In all cases, we have applied a e(re) oc r^^ radial emis- 
sivity and have considered angular emissivities of the form (i) 
/(/ie) oc cons, (solid lines), (ii) f{/J,e) oc /jl'^^ (short dashed lines) 
and (iii) /(/^e) oc (1 + 2.06/ie) (long dashed lines). In the case 
of the Schwarzschild hole, the contribution from the higher order 
images is ^ 20% of the total luminosity, even at high inclinations. 
However, in the case of the extreme Kerr black hole, at inclina- 
tions of ^ 75°, the first order image can dominate the luminosity 
of the system, due to the formation of caustics within the image. 



interfero meter imaging superm assive black holes in nearby 
galaxies (iGendreau et al.ll200iri. 



4 CONCLUSIONS 

Photons orbiting a black hole encounter the strongest pos- 
sible gravitational fields compatible with still escaping to 
the observer. Thus these higher order null geodesies provide 
the best test of Einsteins gravity in the strong field limit. 
We have developed a new strong gravity code capable of de- 
scribing these paths, and calculate them for a geometrically 
thin, optically thick (standard) disc in both Schwarzchild 
and Kerr metrics. These higher order image paths must cross 
the equatorial plane, so are absorbed where this is filled by 
the optically thick disc. As has long been known, the ma- 
jor amplification effects of gravitational lensing are for the 
first order paths from the far side of the underneath of the 
disc viewed at high inclination i.e. photons initially emitted 
downwards on the far side of the black hole, which are bent 
by gravity up above the disc plane. Most of these paths will 
re-intersect the disc unless it has very limited outer rad ial 
extent. While such discs may exist (iRevnolds et a l. 1 1200 A . it 
seems far more likely that these photons would be absorbed 
by material in the equatorial plane. However, there is some 
fraction of the higher order images where the light paths are 
so strongly bent that they re-intersect the equatorial plane 
very close to the photon orbit radius. By definition, this is 
below the minimum stable orbit for particles, so the stan- 
dard disc cannot exist at this point. Instead, for a stress- free 
inner boundary condition, the disc material plunges rapidly 
though this region, so there is much less absorbing mate- 
rial in the equatorial plane. This material may be optically 



thick at high mass accretion rates (|R.evnolds BegelmanI 
119971^ . but this depends on whether the flow in the plung- 
ing region in smooth or clumpy. Hence higher order photons 
which cross the plane below 6rg need not be reabsorbed by 
the disc. We show that these observable higher order photons 
can carry 10% of the flux for a face on disc. Edge on discs 
reduce the expected observable flux as only about half of the 
orbiting photon ring can be seen through the gap below rms , 
while the rest re-intercepts the disc at r ^ rms • 

The situation is less favourable in the Kerr geometry 
as photons now orbit at a range of radii depending on how 
their angular momemtum is aligned with the spin of the 
black hole. Photons going against the spin will orbit at radii 
which are larger than that of the minimum stable orbit of 
the disc, so should be absorbed rather than escaping through 
the inner 'hole' in the disc (provided a > 0.65). Thus for 
Kerr black holes, even though there is a greater fraction of 
photons in the higher order images, we expect that these 
are less observable due to the overlap between the photon 
and particle orbits. Also, even though there is still a gap 
between the last stable particle orbit and the aligned pho- 
ton orbits, this gap is much smaller than in Schwarzchild, 
both in terms of radial coordinate and in terms of impact 
parameter on the sky. For a=0.998, the equilibrium spin for 
thin disc accretion, rms = 1.23r^ while the r^^^ = 1.07, so 
Ar = 0.16. Since the photons orbit in slightly stronger grav- 
ity, their paths are more distorted, so the difference in their 
impact parameter at infinity is slightly less than the dif- 
ference in radial coordinate. Thus the higher order images 
round a rapidly spinning black hole are much more diffi- 
cult to spatially resolve from the direct image of the disc, 
and carry less flux than Schwarzchild as only a fraction re- 
intersect the plane below rms • 

Thus the best chance of observing these higher order 
photon paths from a standard accretion disc is with a face- 
on Schwarzchild black hole. X-ray interferometry could po- 
tentially directly resolve such scales for supermassive black 
holes in nearby galaxies, and such missions ar e being seri- 
ously proposed fMAXTM: fOendreau et al.l2 001^. Even if the 
plunging region is optically thick, direct imaging with res- 
olution ^ 2rg will clearly show the black hole spin from 
the apparent size of the 'hole' in the center of the disc. 
The radial coordinate of the innermost stable orbit in an 
equilibrium spin Kerr metric is a factor of 3 smaller 
than for Schwarzchild, which translates to a change in ap- 
parent size at infinity of the 'hole' diameter from ^ 5r^ 
(a=0.998) to ^ 14r^ (a=0). This is important as recent pa- 
pers have emphasised that the true size of the shadow of 
the event horizon (in effect the impact parameter of the 
orbiting photons, given fairly accurately by our second or- 
der image) is rather similar i n both Schwarzchild and Ker r 
jFalcke. Melia. A^ollEinfl WukuAEm^ iTa.kaha.shilEinl. 
This is true, but different to the size of the 'hole' defined by 
the innermost stable orbit of the disc. For continuous energy 
release down to the event horizon, the 'hole' in the image 
of the accretion flow is set by the true shadow size, but for 
accretion flows where the energy is only released from stable 
particle orbits (the disc, as opposed to the plunging region) 
then the innermost stable orbit sets the size scale. Whether 
the energy release can be contin uous across rms is cur - 
rently a matter of active research (|Krolik Hawlevll2002r i. 
However, there is clear observational evidence for an inner- 
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most stable orbit i n the disc dominanted spectra of galactic 
black hole binaries (lEbisawa et al .11 19931 : iKubota et al.lll999l : 
lOierlihski Pond 1200 A . so it seems likely that the accre- 
tion flow can take the standard thin disc form assumed here. 

The spectral signatures of these photons only become 
apparent for high-inclination systems, where they carry be- 
tween ^ 20% (Schwarzschild) and ^ 60% (extreme Kerr) of 
the total luminosity of the system. However, in a realistic 
system, the disc extends out to many thousands of gravita- 
tional radii from the black hole and so these photons return 
to the disc before reaching the observer. Far from the central 
black hole, the effect of the returning radiation is compara- 
ble to the gravitational potential energy of the material and 
so these photons can play an important role in shaping the 
properties of the disc here. In particular, reprocessing of this 
returning radiation at large distances from the hole will po- 
tentially provide a physical mechanism to flatten (at large 
radii) the r~^^'^ temperature profile of a planar accretion disc 
irradiated by a central source, which could help to resolve 
the conflict between the predicted and observed optical /UV 
colours in Active Galactic Nuclei (iKoratkar Blaeslll999l) . 

The code used to calculate these re sults is an exten- 
sion of that of iBeckwith Pond ()2004ri . The new aspect 
described here is a set of analytic constraints on the pos- 
sible photon trajectories which vastly reduce the scale of 
the calculation. This, combined with the use of the analytic 
integrals (elliptic functions) for the photon paths, has the 
result that the calculation of a pho ton traiectory linking 
two points (the emitter and observer. IViergut3l 19931) can be 
performed by a simple minimisation, capable of being calu- 
lated to almost arbitrary accuracy on a standard desktop 
pc on time scales of a few minutes. This technique allows us 
to include the contribution of orbiting photons in the calcu- 
lations without loss of resolution. The code is also flexible 
enough to be adapted to any accretion flow, and a future pa- 
per will consider the contribution of the higher order images 
in an optically thin, geometrically thick flow. 
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